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C^! Abstract 

D , In models with TeV-scale gravity, ultrahigh energy cosmic rays can generate microscopic black 



holes in the collision with atmospheric and terrestrial nuclei. It has been proposed that stringent 



^ ■ bounds on TeV-scale gravity can be obtained from the absence of neutrino cosmic ray showers 

H : 

5r , mediated by black holes. However, uncertainties in the cross section of black hole formation and, 

most importantly, large uncertainties in the neutrino flux affects these bounds. As long as the 
cosmic neutrino flux remains unknown, the non-observation of neutrino induced showers implies 
less stringent limits than present collider limits. 
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I. INTRODUCTION 

It has been proposed that in models of TeV-scale gravity Q], ultrahigh energy cosmic 
neutrinos colliding with atmospheric nuclei could form black holes (BHs) in the atmosphere 
P]. (For a review on nonperturbative gravitational events at trans-Planckian energies and 
references, see e.g. Ref. The products of BH decay would then be detected as extensive 
air showers (EASs) of ultrahigh energy cosmic rays (UHECRs). 

One of the signatures of BH formation in the atmosphere would be the observation 
of deeply penetrating quasi-horizontal EASs. Semiclassical calculations of BH cross sections 
suggest an interaction length for neutrino-nucleon events in air of the order of 10^ g cm~^, 
about two orders of magnitude lower than the standard model (SM). This is considered 
enough to generate deeply penetrating horizontal air showers with little background from the 
SM. The nonobservation of these deeply penetrating horizontal air showers can set an upper 
bound on the BH cross section . Since the semiclassical BH cross section is inversely 

proportional to the fundamental gravitational scale M^, a lower bound on follows. For 
a number of extra dimensions n > 5, the estimated lower bound on the gravitational scale 
is Md > 1.0-1.4 TeV Q. 

Limits on the fundamental gravitational scale from air showers depend on two main 
assumptions: i) the existence of a cosmogenic neutrino flux; and ii) the accuracy of estimates 
for the cross section for BH formation. The cosmogenic neutrino flux is not known with 
good precision. Many models have been proposed to estimate the flux, which varies by 
more than an order of magnitude (see e.g. 

experiments are quite generous in their range In addition, microscopic BH formation at 
trans-Planckian energies is not understood nor has it been observed in particle collisions. 
Recent work assumes that the cross section for BH formation at parton level is approximately 
the semiclassical black disk (BD) area with Schwarzschild radius rs{M), where M is the 
center of mass (cm.) energy of the collision. However, a rigorous calculation of the BH 
cross section is still an open question. 

The detection of microscopic BH formation in particle collisions would confirm TeV- 
gravity models. However, nonobservation of these events does not necessarily rule out TeV- 
gravity at this stage. A more stringent constraint on M^, cannot be set until the neutrino 
flux and the physical process of BH formation are better understood. Microscopic BHs might 
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not form in the atmosphere even if the gravitational scale is of order of the TeV. Therefore, 
constraints on Mr, from nonobservation of atmospheric BH events are not as stringent as 



what was previous^ 
bounds are based 
gravitational effects. 



considered. The perturbative predictions, on which current collider 
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, seem to be a much stronger basis than those invoking strong 



In this article we revisit and discuss lower bounds on Mo from nonobservation of BH- 
induced EASs in more detail, and find how various uncertainties may reduce these bounds. 

II. CROSS SECTION 

The total cross section for a neutrino-nucleon BH event in {n + 4)-dimensions is obtained 
by summing the neutrino-parton cross section of BH formation aiyi^BH{xs,n, Md) over the 
parton distribution functions (PDFs) qi{x, 



(Tup^bh{s, Xm, n, Md) = X] / dx qi{x, -Q^) a, 



ui^BH{xS,n,MD) , (1) 

where —Q^ is the four-momentum transfer squared, ^/x is the fraction of the nucleon mo- 
mentum carried by the parton, y/sxm = MBH,min, the minimal BH mass where semiclassical 
description is vahd, and the parton cross section (7^i^BH{xs,n, Md) is given by 



(Tui^BHixs, n, Md) = F 



(2 + n) 



4t 2 
'xs \ 



D 



(2) 



where F is a form factor. 

It should be noted that MBH,min is not necessarily equal to the minimum allowed mass of 
the BH, M, and MBH,min ^ M > Md ■ MBH,min and M depend on quantum gravity physics 
and cannot presently be determined. Xm is generally assumed to be a constant parameter of 
order one. For spherically symmetric BHs a justification for this choice is provided by the 
following semiclassical argument j^: For MBH,minl Md ^ 3 and n > 5, the Hawking entropy 
of the BH is larger than 10, and therefore strong gravitational effects can be neglected. 
The semiclassical results based on Eq. are then extrapolated for MBH,min/ Md ^ 3 with 
the assumption that the BH or its Planckian progenitor decays on the brane, whatever the 
quantum theory of gravity may be. As M and Md are not necessarily equal, caution is 
required in extrapolating down the semiclassical results to Mbh, min ~ M^ ~ M^Di whcre 
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the physics is unknown and quantum effects are important [17|. The entropy estimation is 
based on Hawking's semiclassical theory and is not vahd at energies a few times above the 
fundamental scale. For example, it has been shown that the existence of a minimum length 
could dramatically increase the value of M, and thus also increase Xm Q]- Finally, the 
estimate simply refers to static spherically symmetric BHs, and may be drastically affected 
if the geometry of the BH or its Planckian progenitor is different. 

The form factor depends in principle on the energy of the process, on the gravitational 
scale, on the geometry and number of extra dimensions, and on the geometry and physical 
properties of the gravitational object. Some of the physical parameters that can affect 
the form factor are angular momentum, charge, geometry of the trapped surface, quantum 
corrections to classical gravity, unknown effects of super-Planckian particle physics, structure 
and topology of the compactified dimensions. With the lack of further insight, most authors 
simply set F = 1. Yoshino and Nambu (YN) numerically investigated the formation of 
the BH apparent horizon. In the YN approach F is a numerical factor depending on the 
number of extra dimensions ranging from ^ 0.65 {n = 0) to ~ 1.88 {n = 7). YN gives 
a relation between the impact parameter and the mass of the BH which is formed in the 
collision. The result is that the mass of the BH decreases as the impact parameter increases 
up to a maximum value. This behavior affects the computation of the total cross section 
by requiring the lower bound of the integral in Eq. (P) to depend on the impact parameter. 
The total cross section in the YN approach is ^ 



Kp^BH{^^^m,n,MD) = '^zdz j dxqi{x,-Q^)a^i^BH{xs,n,MD) 



(3) 



where z is the impact parameter normalized to its maximum value and x'^ = Xm/y'^{z), y{z) 
being the fraction of cm. energy that is trapped into the BH. The YN approach lowers the 
BD cross section (Fig. P). 

The YN result still relies on a number of assumptions that may affect the final estimate. 
(For recent criticisms, see Ref. I20I.) The incoming partons, for example, are modeled as 
classical neutral point-like particles. Partons carry color and EM charge, and it has been 
shown that the physics of collisions between charged particles is quite different from that 
of uncharged ones 0|. Moreover, it is not clear what constitutes the energy that is not 
trapped inside the horizon. Recent studies seem to_suggest that gravitational emission can 
account only for a part of the missing BH energy ^] . This could signal that the physics of 
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FIG. 1: BH cross section {MBH,min = Md = 1-0 TeV) for BD approximation (left upper curves) 
and for YN formalism (right lower curves). For each group, n = 3 ... 7 from below. Q = and 
and CTEQ6 PDFs is used. 

trans-Planckian collisions is more complex than the simple semiclassical picture. 

Another source of uncertainty in Eq. and Eq. comes from the ambiguity in the 
definition of the momentum transfer for a BH event j23|. The latter is usually chosen to 
equal the BH mass or the inverse of the Schwarzschild radius. However, there are no definite 
arguments to prefer either one or to exclude alternative choices. The uncertainty due to the 
ambiguity in the definition of the momentum transfer is evaluated as ~ 10 — 20% jsl. 

Finally, a minor but additional source of uncertainty in the total cross section is due 
to the PDFs. Different sets of PDFs are defined in the literature. The PDFs are not 
known for momentum transfer higher than a given value (see, e.g. Ref. which is lower 

than the momentum transfer expected in BH formation. Equation (Q) and Eq. are 
calculated by imposing a cut-off at this energy. They also suffer from uncertainties at any 
momentum transfer that can contribute to the reduction in the total cross section. The 
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minimum uncertainty on the BH total cross section due to the PDFs can be estimated for 
a given distribution. The CTEQ6 distribution gives an uncertainty of ~ 3 — 4%, a value 
that does not include the uncertainty due to the cutoff on the momentum transfer nor the 
uncertainty introduced by the use of different sets of PDFs. The MRST distribution [25] for 
M^) = 1.0 TeV and MBH,min/MD = 3 gives results about ~ 10 — 15% lower than the CTEQ 
distribution. 



III. EFFECT OF CROSS SECTION UNCERTAINTIES ON Md BOUNDS 

The BD approximation may be different than the actual cross section due to these uncer- 
tainties. To give a concrete example of how the determination of the cross section affects the 
Md bounds, let us look at a case where the parton cross section is arbitrarily reduced by a 
constant factor of order one. The cause of this reduction could be any one or a combination 
of the uncertainties previously discussed. 

Consider high energy primary neutrinos with energy E^, = 10^ — 10^ TeV, MBH,min = 3M£) 
and n = 5. For the sake of simplicity, we temporarily neglect the YN results and consider 
Eq. (^. This is sufficient for a rough estimation. The total cross section (0) for Md = 1 
TeV and F = 1 lies within the shaded band of the total cross section for M^, = 0.72 TeV 
and F = 1 ± 2/3. This is shown graphically in Fig. |21 Therefore, what is interpreted as the 
Md = 1.0 TeV bound in the naive BD approximation can actually be a less stringent bound if 
the parton cross section is two-thirds smaller: Md = 0.72 TeV, the lower bound from collider 



experiments derived from virtual graviton exchange [Ij, The situation is made even 
more complicated by the presence of the unknown parameter MBH,min- Increasing MBH,m.in 
corresponds to decreasing the BH cross section. The lower bound on the eleven-dimensional 
Md for MBH,min/ Md = 5 is approximately half the bound for MBH,min/ Md = 1- 



IV. NUMBER OF EVENTS AND NEUTRINO FLUX 

The uncertainties in the BH cross section listed above apply to BH formation by UHECR 
events as well as by particle colliders. In cosmic ray events, however, the unknown neutrino 
flux adds further uncertainties. In the YN approach, the number of neutrino-nucleon BH 
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Energy (TeV) 

FIG. 2: The cross section for M^, = 1.0 TeV (thick sohd hne) and for Mf, = 0.72 TeV (dashed 
hne). Shaded band denotes F = 1 it 2/3 for Md = 0.72 TeV. The cross section of = 1 TeV 
with F = 1 lies within this band. CTEQ6 is used. 



events detected by a cosmic ray detector in time T is 



dE^ / 2zdz / dxqi{x,-Q'^)ayi^BH{xs,n,MD) —- 

•-'0 J X" V 



A{yE,) , (4) 



where A{yEy) is the experiment acceptance for an air shower energy yE^, Na is Avogadro's 
number, and (i$ / dEi, is the source flux of neutrinos. 

The cosmogenic neutrino flux is considered to be the most reliable source of neutrinos. 
In this model j26| neutrinos are produced from ultrahigh energy protons interacting with 
the ubiquitous cosmic microwave background. However, this is not fully guaranteed as the 
existence of the cosmogenic neutrino flux relies on the assumption that cosmic rays with 
energies above 10® TeV are extragalactic protons. Neither the source nor composition of 
cosmic rays above 10® TeV are known (see, e.g., Ref. 0,Q) 

If these are heavy nuclei or 

photons, or Lorentz invariance is violated [2^ , there may be no cosmogenic neutrino flux at 
ultrahigh energies. Even if UHECRs are protons from extragalactic sources, cosmological 
evolution, spatial distribution, abundance, and injection spectrum of UHECR sources can 



change the cosmogenic neutrino flux by an order of magnitude. In order to derive a conser- 
vative lower bound on Md, the lower end of plausible cosmogenic neutrino fluxes should be 
used in Eq. (H)). 

Here we consider a model that gives a relatively low neutrino flux ii i ag reement with 



observations ^. The flux is calculated following the procedure of Ref. The source 

spectrum is proportional to x exp{—E/Ec), where Ec = 10^'^ TeV is the cutoff energy, 

and normalized to the cosmic ray luminosity at 10^ TeV. There is no cosmological evolution 
and the redshift integration is from z = 0.05 — 8.00. The parameters used are consistent 
with observations: A spectral index of 2.6 is the best flt to the highest energy data and 
the cutoff energy is limited by the highest energy cosmic ray event observed to date. Since 
the ultrahigh energy proton sources are unknown, their evolution cannot be determined. No 
evolution is a reasonable assumption for the lower end of the neutrino flux that should be 
considered. Our model flux is in good agreement with the lower bound obtained from cosmic 
ray data analysis . Figure |H1 compares this flux to two other fluxes, by Protheroe and 
Johnson (PJ) [l^ and Engel, Seckel, and Stanev (ESS) 11]. Our flux is lower by at least 
an order of magnitude than the PJ flux on the most relevant range of neutrino energies. 
The sensitivity of the detector for horizontal air showers also affects the computation of 



the lower bounds on M^- For example, the Pierre Auger Observatory [31[ uses 1.2 m high 
water Cherenkov detectors whereas AGASA 3^ uses 5 cm thick scintillators. From purely 
geometrical arguements, AGASA's detection capability rapidly goes down at large zenith 
angle. On the other hand, AGASA has a lower detector trigger threshold which gives higher 
sensitivity at lower energies. Scaling the aperture by the experiment size is an approximation 
that should be further improved. 



V. REALISTIC BOUNDS ON Md 

The effect of the cosmogenic neutrino flux on M]j lower bound is best illustrated by 
comparing the conservative flux of Sect. IV and the PJ flux. We performed a systematic 
analysis of the number of events using the above cosmogenic neutrino fluxes for different 
experiments (AGASA [33], HiRes y, RICE Q) 

and different choices of parameters in 

^ We thank Todor Stanev for providing the data. 
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FIG. 3: Cosmogenic neutrino fluxes for the t'^ family. PJ (dashed line), ESS (dotted line), and our 
model (solid line) fluxes from the top. 

the total cross sections from Eq. (HI) and (01). For simplicity we used the apertures eiven 
by Refs. |6|, 17|, |8|| . Though RICE is an experiment looking into "ice showers" rather than 
EASs, the technique used in EAS experiments is applicable. The use of the PJ flux gives a 
lower bound of M^i > 1.0 — 1.4 TeV, which is comparable to and sometimes more stringent 
than collider bounds. The most optimistic case with our model flux is when Q = Mbh 
and MBH,min = The lower bounds on the fundamental scale range from = 0.3 

{MBH,min = Mb) to Mn < 0.2 TeV {MBH,min = 5Mb)- These limits are much lower than 
the collider limits ^ jl^-list . 

We also considered changing the cosmological evolution of the conservative flux to be 



^ The computations were performed using the CTEQ6 PDFs and YN approximation. The use of MRST 
PDFs give results identical within an uncertainty of 0.05 TeV and the BD approximation increases Mu 
by 0.05-0.15 TeV. 
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that of Ref. 36|, i.e. 



(1 + 1.9)3 

(1 + 1.9)3exp[(2.7-z)/2.7] 



{z < 1.9) 

(1.9 <z< 2.7) 

(2.7 < z) 



(5) 



This increases the flux by about an order of magnitude. It is lower than the ESS flux for 
^ 10^ TeV, and is considerably less by a couple of factors than the PJ flux in all energy 
ranges. The lower bounds increase by about 0.4 TeV, but are still below the collider bounds. 

VI. CONCLUSIONS 

We have reviewed the method of constraining the fundamental Planck scale from nonob- 
servation of EASs and discussed the uncertainties on the Md lower bound coming from the 
BH cross section and neutrino flux. The major source of uncertainties in the cross section is 
due to the lack of a definite theory of trans-Planckian scattering. Trans-Planckian gravity 
is not known at present to formulate a reliable model for BH formation. The BD approx- 
imation and its variants seem to provide a reasonable model for the BH formation cross 
section at parton level. However, the naive BD cross section relies on a number of crude 
assumptions and suffers from uncertainties that cannot be estimated. Neglected quantum 
gravity effects are expected to be relevant. 

A conservative estimation of Md lower bounds is obtained by taking into account these 
caveats and considering the range of cosmogenic neutrino fluxes which are compatible with 
observations. The conclusion is that in the absence of independent determination of the 
cosmogenic neutrino flux and of reliable theoretical computation of the BH cross section, 
nonobservation of BH-induced deeply penetrating EASs gives lower bounds on the funda- 
mental gravitational scale less stringent than present collider bounds. 
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Note added in Proof 

Recent work by Hossenfelder [^^l calculates the BH cross section suppression from min- 
imal length at LHC. It is found that the ratio between the total cross section with and 
without minimal length effects is approximately 0.19 for the expected LHC energy, and in- 
creasing at higher energies. This result strenghtens our arguments of section II, namely that 
the presence of a minimum length may substantially reduce the BD cross section. 

In addition, progress in neutrino flux estimates with cosmic ray primaries other than 
protons have shown how uncertain the neutrino flux can be (see |hep-ph/0409"3T6| and 
hep-ph/0407"618l ). 
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